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We have employed first principles calculations to study the electronic structure and magnetic 
properties of the low- dimensional phosphates, Ba2Cu(P04)2 and Sr2Cu(P04)2- Using the self- 
consistent tight binding linearized muffin-tin orbital (TB-LMTO) method and the N th order muffin 
tin orbital (NMTO) method we have calculated the various intrachain as well as the interchain 
hopping parameters between the magnetic ions (Cu 2+ ) for both the compounds. We find that the 
nearest-neighbor intrachain hopping (t) is the dominant interaction suggesting the compounds to 
be indeed one-dimensional. Our analysis of the band dispersion, orbital projected band structures 
and the hopping parameters confirms that the Cu 2+ - Cu 2+ super-super exchange interaction takes 
place along the crystallographic 6-direction mediated by O-P-O. We have also analyzed in details 
the origin of short range exchange interaction for these systems. Our ab-initio estimate of the ratio 
of the exchange interaction of Sr2Cu(P04)2 to that of Ba2Cu(P04)2 compares excellently with 
available experimental results. 



PACS numbers: 71.20.-b, 71.70.Gm, 75.10.Pq 
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I. INTRODUCTION 

Low-dimensional quantum spin systems with chain, 
ladder, or planar geometries have attracted much atten- 
tion due to their unconventional magnetic properties. 1 
Much effort has been devoted particularly over the 
last decade to understand the behavior of quasi-one- 
dimensional spin systems. These systems exhibit a rich 
variety of phases due to the enhanced quantum fluctu- 
ations in reduced dimensionality. It is well known that 
the half-integer uniform antiferromagnetic chain, though 
having a singlet ground state, is gapless and quantum 
critical with power law spin correlations. A deviation 
from the uniform chain due to, for instance, an alter- 
nating exchange leads to a gapped (spin-gap A)ground 
state. Frustration arising in one dimensional spin chains 
due to next nearest neighbor interactions can also lead to 
a gap in the spin excitation spectrum. A well known case 
is that of the Majumdar-Ghosh(MG) chain^ where the 
Hamiltonian is antiferromagnetic in nature, and includes 
both nearest neighbor J\ and next nearest neighbor J 2 
exchange interactions. This model is exactly solvable for 
the ratio ^ = a = 0.5. The MG-chain supports a spin 
gap for a > a cr — 0.2411. For < a < a cr , a gapless 
phase is obtained. In contrast, integer-spin uniform AF 
chains are always gapped as conjectured by Haldane 3 as 
also half-integer even-leg ladders^ while the half-integer 
odd- leg ladders are gapless. Recently, the discovery 
of spin-Peierls transition in the one-dimensional (ID) 
Heisenberg antiferromagnet (HAF) CuGeOg^ renewed 
interest in 1-D transition metal oxides with spin S — 1/2 
ions such as Cu, which are responsible for both the mag- 
netic properties as well as the Jahn- Teller distortion of 
the lattice. 

Low-dimensional quantum spin systems are important 
also due to the belief that understanding these systems 
may be crucial to understand the physics of high-T c 



cuprates^ In addition, the field of low-dimensional quan- 
tum magnetism provides a fertile ground for rigorous the- 
ory. Powerful techniques like the Bethe Ansatz (BA) and 
bosonization are available to study ground and excited 
state properties of these systems. Models of interact- 
ing spin systems are known for which the ground state 
and in some cases their low-lying excitation spectra are 
known exactly. The knowledge gained from these mod- 
els provides impetus to look for real materials so that 
experimental confirmation of the theoretical predictions 
can be made. 

Quasi-ID magnets, in which the magnetic interaction 
in one direction dominates, with much weaker interac- 
tions in other directions, exhibit short-range order over 
a wide temperature range. A uniform Heisenberg half- 
integer-spin chain as discussed above has a gapless spin 
excitation spectrum. The ground state does not have 
long range order (LRO) because of strong quantum fluc- 
tuations. In real systems however, inter-chain interac- 
tions can lead to a deviation from ideal behavior and 
even magnetic LRO. Whereas a large number of quasi- 
1D systems have been studied, there are only a few 
cases in which the interchain interactions are very weak 
and true ID behaviour can be observed down to low- 
temperatures. One example is Sr2Cu03 where the ex- 
change coupling constant J/ks is about 2200±200 K and 
the ordering temperature Tjv = 5.4 K, giving the ratio 
k B T N /J w 0.25%1 Another example is copper benzoate 
Cu(C 6 H 5 C00) 2 -3H 2 with J/k B = 17.2 K4 

Very recently, the susceptibility measurements by Belik 
et al£ and susceptibility and nuclear magnetic resonance 
(NMR) measurements by Nath et al^ have suggested 
that Ba 2 Cu(P0 4 )2 (J/k B = 151 K) and Sr 2 Cu(P0 4 ) 2 
(J/k B = 165 K) are excellent S — 1/2 linear chain 
Heisenberg antiferromagnets. The temperature depen- 
dence of the NMR shift K(T) is well described by the 
5=1/2 Heisenberg antiferromagnetic chain model. 11 



The J/ks in these compounds is about one order of 
magnitude smaller than in Sr 2 Cu03 and about one or- 
der of magnitude larger than in copper benzoate. In 
Sr2Cu(P04)2, Bclik et a/£ suggested from their mea- 
surements that LRO sets in at TV = 0.085 K, which 
implies that IzbTn/J = 0.06%, making Sr2Cu(PC>4)2 the 
best realization of a ID S — \ HAF system studied to 
date. This provides a rare opportunity for a compari- 
son of experimental results with theoretical models and 
improving our understanding on these systems. 

The unique ID magnetic behavior exhibited by these 
compounds may be traced back to the interplay of the 
geometry (crystal structure) and quantum chemistry and 
therefore requires a thorough understanding of the elec- 
tronic structure of these systems. In this paper, we 
shall study the electronic structure of Ba2Cu(P04)2 and 
Sr2Cu(PC>4)2 in some detail. The characteristic feature 
of these compounds is a chain of isolated CUO4 square 
plaquettes along the crystallographic b direction. The 
goal of the present work is to understand and evaluate 
the various exchange paths leading to ID magnetism by 
calculating the various hopping integrals. The latter will 
also provide the dominant exchange integrals entering the 
magnetic Hamiltonian for the spin-half Heisenberg chain 
given by 
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FIG. 1: (Color online) (a) Schematic diagram of [Cu(PC>4)2]oo 
linear chains propagating along the b direction for 
(Ba/Sr)2Cu(P04)2- (b) Various hoppings between the Cu 
atoms, (c) Schematic energy level diagram for a Cu d ion in 
square planar environment 



where Jy is the magnetic exchange interaction and S 
denotes the spin. The various hopping matrix elements 
between the magnetic ions (Cu for these compounds) are 
strongly dependent on the way the CUO4 squares are 
assembled. The edge sharing CUO4 squares and corner 
sharing CUO4 squares lead to quite different values for 
the nearest-neighbor hopping (t) and the exchange inter- 
action (J). In the latter case, the 180° Cu-O-Cu bond 
with a common O 2p orbital gives rise to large magni- 
tude of J as seen for the Cu02 planes in the high-T c 
cuprates. This contrasts with the 90° Cu-O-Cu bond 
when two CUO4 squares share an edge, where the hopping 
matrix element t and therefore the exchange interaction 
are expected to be small in magnitude. In the present 
work, we shall explicitly obtain the various effective Cu- 
Cu hopping parameters using N-th order muffin-tin or- 
bital method (NMTO) and clarify the role of electronic 
structure on the magnetic properties of these systems. 

The remainder of the paper is organized as follows: 
In Sec. II, we describe the crystal structure of both 
Ba2Cu(P04)2 and Sr 2 Cu(P04)2 followed by the compu- 
tational details. Section III is devoted to the detailed 
discussion of the electronic structure using tight-binding 
linearized muffin-tin orbital method (TB-LMTO) within 
atomic sphere approximation (ASA) and the NMTO 
method followed by a comparative study of our calcu- 
lations with available experimental results. Finally, in 
Sec. IV, we present our conclusions. 



II. CRYSTAL STRUCTURE AND 
COMPUTATIONAL DETAILS 

Ba2Cu(P04)2 and Sr2Cu(P04)2 are two isostructural 
compounds in the family of S = 1/2 ID HAF systems. 12 
They both crystallize in the monoclinic structure with 
space group C 2 / m (No. 12). The lattice parameters 
are a = 12.16 A, b = 5.13 A, c = 6.88 A, and = 
105.42° and a = 11.51 A, b = 5.07 A, c = 6.57 A, and 
13 = 106.35° for Ba 2 Cu(P0 4 )2 and Sr 2 Cu(P0 4 )2, respec- 
tively. The characteristic features of the structure as il- 
lustrated in Fig. 1 are isolated, square, CUO4 plaquettes 
sharing their edges with two similar kind of PO4 tetra- 
hedra. This edge sharing via PO4 tetrahedra takes place 
along one crystallographic direction (b direction) forming 
isolated [Cu(P04)2]oo chains along the crystallographic 
b direction. The chains, are however, staggered with re- 
spect to each other along the crystallographic a direction. 
The Ba/Sr cations reside between these parallel chains. 
In Fig. 1(b), we have indicated the hoppings between 
various Cu atoms. In Fig. 1(c), we have displayed the 
schematic energy level diagram due to crystal field split- 
ting of the Cu d ion in the square planar coordination 
provided by the oxygens. We find that the highest en- 
ergy level is Cu d x 2_ y 2 and is well separated from the 
rest of Cu d levels. 

Our analysis of the electronic structure is carried out 
in the framework of the (TB-LMTO) methodic in 
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the (ASA) within local density approximation (LDA) to 
the density functional theory. The space filling in the 
ASA is achieved by inserting empty spheres and by in- 
flating the atom-centered nonoverlapping spheres. The 
atomic radii are chosen in such a way that (i) the charge 
on the empty spheres is negligible and (ii) the overlap 
of the interstitial with the interstitial, atomic with the 
atomic, and interstitial with the atomic spheres remains 
within the permissible limit of the ASA. The basis set 
for the self-consistent electronic structure calculation for 
Ba 2 Cu(P04) 2 includes Cu(s, p, d), Ba(s, d and /), 0(s, 
p), and P(s, p), and for Sr 2 Cu(PC>4)2 the basis set in- 
cludes Cu(s, p, d), Sr(s, d), 0(s, p), and P(s, p), the rest 
being downfolded. The (8, 8, 8) k mesh has been used 
for self-consistency. All the fc-space integrations were 
performed using the tetrahedron method. 1 - In order to 
extract the various hopping integrals we have employed 
NMTO based downfolding method The downfold- 
ing method consists of deriving a few-orbital effective 
Hamiltonian from the full LDA Hamiltonian by integrat- 
ing out high-energy degrees of freedom. 



III. RESULTS AND DISCUSSION 




-6 -4-2 2 
ENERGY (eV) 



4 
3.5 

3 
2,5 

2 
15 

1 

0.5 




: ( C ) 










■ 01+Q2 


i 1 1 


. 03 ; 








1 -0.5 





0.5 





(d) 


s ! | 



















4 
3.5 

3 
2,5 

2 

15 

1 

0.5 





-4-2 2 
ENERGY (eV) 



: (f) 


1 1 

|Ef\ 




■Ol+O^ 


. —r \ 
\ i 


^.03 ; 






1 -0.5 





0.5 



ENERGY (eV) 



ENERGY (eV) 



The all-orbital band structure and the total density 
of states (DOS) and partial DOS for Ba 2 Cu(P0 4 ) 2 and 
Sr 2 Cu(P04) 2 are displayed in Fig. 2. The bands are 
plotted along the various high symmetry points of the 
Brillouin zone corresponding to the monoclinic lattice. 
All the energies are measured with respect to the Fermi 
level of the compound. The characteristic feature of the 
non-spin-polarized band structure displayed in Figs. 2(a) 
and 2(d) is an isolated half- filled band at the Fermi level. 
This band is predominantly derived from the antibonding 
linear combination of Cu d x 2_ y 2 and O p a states residing 
in the same square plaquette. The conduction band is 
separated from the other Cu d character dominated va- 
lence bands by a small gap. This is consistent with the 
schematic energy level diagram displayed in Fig. 1(c). 
Further below these Cu d bands, are the bands with pre- 
dominantly oxygen 2p character. The conduction band 
disperses only along the chain direction with hardly any 
dispersion perpendicular to the chains. This is reflected 
in the partial DOS shown in the inset of Figs. 2b and 
2e, which has a characteristic feature of ID tight-binding 
dispersion. 

This isolated Cu d x 2_ y s conduction band is respon- 
sible for the low-energy physics of this material. This 
band hardly couples with any other states except for the 
oxygens (03) residing in the square plaquette as can be 
seen from the plot of the oxygen partial density of states 
displayed in Figs. 2(c) and 2(f). Since the calculations 
are carried out in the framework of LDA, the system 
is not insulating. We have checked that the insulating 
behavior can be recovered by including the Coulomb in- 
teraction in the framework of LDA + U. In the follow- 
ing we shall employ NMTO downfolding method to map 



FIG. 2: Band structure, density of states, and partial DOS of 
[(a)-(c)] Ba 2 Cu(P0 4 )2 and [(d)-(f)] Sr 2 Cu(P0 4 )2. 



our LDA results to a low-energy orthogonal tight-binding 
Hamiltonian. This Hamiltonian will serve as the single 
electron part of the correlated Hubbard model relevant 
for this system, which in turn can be mapped to an ex- 
tended Heisenberg model [Eqn. (1)] with the exchange 
couplings related to the LDA hoppings by the relation 



4H 



Jii ~U, 



eff 



(2) 



where U e ff is the screened Coulomb interaction. 

The NMTO downfolding method provides an ab initio 
scheme to construct a low-energy tight-binding Hamilto- 
nian starting from the full all-orbital LDA calculations. 
This method relies on the energy-selective downfolding 
process to integrate out high-energy degrees of freedom 
resulting in a low-energy Hamiltonian defined in the ba- 
sis of effective orbitals. These effective orbitals, by con- 
struction, are tailored to contain in their tails the states 
which are integrated out with weights being proportional 
to their admixture with the orbitals that are retained in 
the basis. If the low-energy set of bands are separated 
from all other bands, the orthonormalized NMTO set 
converges to a set of Wannier functions. Hence, with 
the NMTO downfolding method, we can generate local- 
ized Wannier functions directly without recourse to Bloch 
functions. These effective orbitals (NMTOs) provide a 
direct visualization of the various interaction paths. The 
Fourier transform of this few orbital downfolded Hamilto- 
nian provides various hopping integrals between these 
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effective orbitals, and the corresponding tight-binding 
Hamiltonian can be written as 



(ct Cj -+H.c), 



(3) 



where i and j denote a pair of orbitals. These hop- 
ping integrals thus form the first principles set of param- 
eters obtained without any fitting procedure containing 
the signature of the pathways involved in the hopping 
process. 

For the present compounds with one predominantly Cu 
d x 2_ y 2 band so well separated from the rest, one can have 
a tight-binding model with a single Cu d x 2_ y 2 orbital per 
Cu to be a good approximation to the full band structure 
close to the Fermi level. In Fig. 3, we have displayed 
the band structure of Ba2Cu(P04)2 and Sr2Cu(P04)2 
decorated with Cu d x 2_ y 2 character together with the 
downfolded band structure where only Cu d x 2_ y 2 or- 
bital is retained in the basis and the rest are downfolded. 
The agreement between the two is remarkable, suggesting 
that our one band Hamiltonian is indeed physically rea- 
sonable and is suitable to capture the low-energy physics 
of the system. The Fourier transform of the downfolded 
Hamiltonian H(k) — > H(R) gives the effective hopping 
parameters for the physical Hamiltonian. In Table I, we 
have displayed the various dominant effective hopping 
integrals Uj (having magnitude > 1 meV) between the 
Cu 2+ ions at sites i and j. The notations for various 
hoppings for a pair of Cu site (i and j) are indicated in 
Fig. 1(b), where, e.g., for a Cu ion at site i, = t\ if the 
other Cu ion at site j is at the nearest neighbor-position 
along the chain and similarly for all the other hoppings. 

From Table I, we gather that the nearest-neighbor hop- 
ping (ti) is dominant for both the compounds. The 
other hoppings (t2) along the chain (intrachain) as well 
as between the chains (interchain) (tf, tf) and between 
chains in the perpendicular (c) direction (t±) are an or- 
der of magnitude smaller than the dominant nearest- 
neighbor hopping. The smaller magnitude of the hop- 
ping for Ba2Cu(P04)2 may be attributed to the rela- 
tively large lattice constant for the compound. Recently, 
these hoppings were also calculated by adopting a fitting 
procedure^ Although the trends among the hoppings 
obtained by fitting the conduction band agree with our 
ab initio results, the magnitude as well as the range of 
the hoppings are exaggerated in the fitting scheme. 

In Fig. 4, we have displayed the Cu d x 2_ y 2 Wannier 
function. We note that for both the compounds, the 



Cu d T 



Wannier function is primarily localized in one 



plaquette and it hardly couples with the neighboring pla- 
quettes. The tails of the Cu d x 2_ y 2 orbital are shaped 
according to O p x /p y orbitals such that the Cu d x 2_ y 2 or- 
bital forms strong pdcr antibonds with the O p x /p y tails. 
We also note that the Cu d x 2_ y 2 Wannier function for 
Ba2Cu(P04)2 is more localized, explaining the smaller 
magnitude of the hoppings for this system. 
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projected as well as downfolded band in 



[(a) and (b)] Ba 2 Cu(P0 4 )2 and [(c) and (d)] Sr 2 Cu(P0 4 )2. 
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FIG. 4: (Color online) Effective Cu~d x 2 _ y 2 Wannier function 
plot for (a) Ba 2 Cu(P0 4 )2 and (b) Sr 2 Cu(P0 4 )2. The d x 2_ y2 
orbital is denned with the choice of local co-ordinate system 
as discussed in the text. The spheres represent the ions as 
indicated in the figure. 



The hopping parameters listed in Table I can be used 
to estimate various exchange interactions using Eq.(2). 
In the absence of a satisfactory way of direct computa- 
tion of exchange integrals, such an approximate method 
is a good starting point for estimating the exchange cou- 
plings as well as the relative strengths of the various 
exchange integrals. The ratios of the various exchange 
interactions for each compound are listed in Table II. 
As expected, the nearest-neighbor exchange interaction 
is dominant for both the compounds with practically neg- 
ligible interchain and further-neighbor intrachain inter- 
actions, making these compounds ideal ID HAFs span- 
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ning a wide temperature range. The agreement between 
the ratio of the nearest-neighbor exchange interaction in 
Ba 2 Cu(P0 4 ) 2 to that in Sr 2 Cu(P0 4 ) 2 , 7§f, obtained in 
our calculation and that obtained from experiment is re- 
markable. In the following, we shall try to understand 
the origin of the short-range hopping and therefore the 
exchange interactions in these systems. From the Cu 
d x 2_ y 2 Wannier function (Fig. 4), we gather that the 
Cu-Cu hopping primarily proceeds via the oxygens. As a 
result, as argued by Koo et. aL— the strength of the Cu- 
O...Cu spin-exchange is primarily governed by the 0...0 
distance and Cu-O..Cu angles rather than the Cu..Cu dis- 
tances. The exchange interaction therefore becomes neg- 
ligible when the relative separation between the oxygens 
mediating the interaction is large compared to the van 
der Waals distance. For the present compounds, the weak 
hopping integrals and therefore the exchange interactions 
between the Cu ions may be attributed to their pecu- 
liar geometry. The CUO4 plaquettes in the present com- 
pounds neither share edges nor corners and are connected 
by PO4 tetrahedra resulting in a super-super exchange 
process to mediate the exchange interactions. This ar- 
rangement results in a large relative separation (greater 
than the van der Waals distance) between the oxygens ex- 
cept for the oxygens in the neighboring plaquettes along 
the chain, resulting in an appreciable nearest-neighbor 
interaction while other intrachain as well as interchain 
interactions are negligible. 

As mentioned earlier, even tiny interchain interactions 
can cause a deviation from the ideal ID behavior. In- 
deed, in the 31 P NMR shift (which is proportional to 
the spin susceptibility) measurements on Ba 2 Cu(PC>4) 2 
and Sr 2 Cu(P04) 2 ji£ it was found that while the data are 
well described by the ID HAF model in a large temper- 
ature range, a sharp deviation was seen below ksT/J w 
0.03 (still above the LRO temperature found by Belik et 
al£). A small interchain interaction can therefore be 
of significance at very low temperatures. In order to 
check that small interchain hoppings are indeed impor- 
tant, we calculated the zero-temperature spin-spin corre- 
lation function for a pair of staggered S = \ Heisenberg 
chains [see Fig. 1(b)] (24 sites, with periodic boundary 
conditions) using the exact diagonalization method. We 
confirmed that in the absence of interchain coupling, the 
spin-spin correlation function decays with a characteris- 
tic power-law behavior. However, inclusion of a very tiny 
interchain coupling ^ w 0.02 causes the spin-spin cor- 
relation function to deviate from the power-law behavior 
and acquire an exponential decay. This might then be 
responsible for the sharp decrease in the 31 P NMR shift 
seen at low-temperatures by Nath et alJ£ 

IV. SUMMARY AND CONCLUSIONS 

We have employed TB-LMTO ASA and NMTO 
method to analyse in detail the electronic structure 



of two isostructural phosphates, Ba 2 Cu(P04) 2 and 



TABLE I: Hopping integrals obtained from our ab initio anal- 
ysis for Ba 2 Cu(P0 4 )2 & Sr 2 Cu(P0 4 ) 2 . 





ti (meV) 


t 2 (meV) 


1 ic 1 ic 1 
tl t 2 tj_ 


Ba 2 Cu(P0 4 


) 2 96 


4 


5 12 


Sr 2 Cu(P0 4 ) 


i 2 103 


9 


14 1 


TABLE II: 


Ratios of the 


exchange 


coupling constants 



and comparison with experiment for Ba 2 Cu(P0 4 ) 2 & 
Sr 2 Cu(P0 4 ) 2 . 





J2/J1 


Jlic/ Jl 


JBa/Jsr a 


Jbo./ Jst 


Ba 2 Cu(P0 4 ) 2 


0.0017 


0.0027 


0.869 


0.915 


Sr 2 Cu(P0 4 ) 2 


0.0076 


0.0184 







"Ratio of J /K B of Ba 2 Cu(P0 4 ) 2 to that of Sr 2 Cu(P0 4 ) 2 ob- 
tained from theory. 

'Ratio of J /K B of Ba 2 Cu(P0 4 ) 2 to that of Sr 2 Cu(P0 4 ) 2 ob- 
tained from experiment — 



Sr 2 Cu(P0 4 ) 2 . The characteristic feature of the LDA 
electronic structure of these compounds is an isolated 
half-filled band at the Fermi level derived from the an- 
tibonding linear combination of Cu d x 2_ y 2 and oxygen 
p a states residing on Cu0 4 square plaquettes. The var- 
ious hopping integrals obtained using NMTO downfold- 
ing method suggest that these compounds are indeed one 
dimensional with dominant nearest-neighbor hopping in- 
tegrals along the chains. Our estimate of the exchange 
interaction, in particular the ratio of the exchange in- 
teraction of Sr 2 Cu(P04) 2 to that of Ba 2 Cu(P04) 2 , com- 
pares excellently with the experimental results. The plot 
of the Cu d x 2_ y 2 Wannier functions suggest that the ex- 
change interactions between the magnetic Cu 2+ ions is 
primarily mediated by the oxygens. Our calculations 
confirm that the unique geometry of the CUO4 plaque- 
ttes, which are neither face sharing nor corner sharing, 
is responsible for the unique ID magnetic properties for 
these systems. Finally, we argue that although the inter- 
chain couplings are small, they may still be important for 
the low-temperature properties for this system and may 
explain the sharp deviation from the ID HAF model at 
low-temperature seen in the 31 P NMR shift for these sys- 
tems. 
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